Quantum phases with differing computational power 
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The observation that concepts from quantum information has generated many alternative in- 
dicators of quantum phase transitions hints that quantum phase transitions possess operational 
significance with respect to the processing of quantum information. Yet, studies on whether such 
transitions lead to quantum phases that differ in their capacity to process information remain lim- 
ited. Here We show that there exist quantum phase transitions that cause a distinct qualitative 
change in our ability to simulate certain quantum systems under perturbation of an external field 
by local operations and classical communication. In particular, we show that in certain quantum 
phases of the XY model, adiabatic perturbations of the external magnetic field can be simulated by 
local spin operations, whereas the resulting effect within other phases results in coherent non-local 
interactions. We discuss the potential implications to adiabatic quantum computation, where a 
computational advantage exists only when adiabatic perturbation results in coherent multi-body 
interactions. 

PACS numbers: 



INTRODUCTION 

The study of quantum phase transitions has greatly 
benefited from developments in quantum information 
theory [2|. We know, for example, that the extremum 
points of entanglement and other related correlations co- 
incide with phase transition points 0-0|, and that differ- 
ent phases may feature differing fidelity between neigh- 
boring states [sMlSj. These observations have helped pi- 
oneer many alternative indicators of phase transitions, 
allowing the tools of quantum information science to 
be harnessed in the analysis of quantum many body 
systems[l, 2]. The reverse, however, remains understud- 
ied. If the concepts of quantum information processing 
have such relevance to the study of quantum phase tran- 
sitions, one would expect that systems undergoing quan- 
tum phase transition would also exhibit different oper- 
ational properties from the perspective of information 
processing. Yet, there remains little insight on how such 
relations are apply to quantum information and compu- 
tation. 

In this paper, we demonstrate using the XY model that 
different quantum phases have distinct operational sig- 
nificance with respect to quantum information process- 
ing. We reveal that the differential local convertibility 
of ground states undergoes distinct qualitative change at 
points of phase transition. By differential local convert- 
ibility of ground states, we refer to the following (Fig. 1): 
A given physical system with an adjustable external pa- 
rameter g is partitioned into two parties, Alice and Bob. 



Each party is limited to local operations on their subsys- 
tems (which we call A and B) and classical inter-party 
communication, i.e., local operations and classical com- 
munication (LOCC). The question is: can the effect on 
the ground state caused by adiabatic perturbation of g be 
achieved through LOCC by Alice and Bob? Differential 
local convertibility of ground states is significant. Should 
LOCC operations between Alice and Bob be capable of 
simulating a particular physical process, then such a pro- 
cess is of limited computational power, i.e., it is incapable 
of generating any quantum coherence between A and B. 

We make use of the most powerful notion of differential 
local convertibility, that of LOCC operations together 
with assisted entanglement [13, [lH • Given some infinites- 
imal A, let \G{g))AB and \G{g + A))ab be the ground 
states of the given system when the external parameter is 
set to be g and ^ + A respectively. The necessary and suf- 
ficient conditions for local conversion between \G{g)) ab 
and \G{g + A))ab is given by Sa{g) > Sa{g + A) for all 
a, where 

Sa{g) = j^iog^[Trp%{g)] = -^log^ [^Af] (1) 

1 — a 1 — a 

1=1 

is the Renyi entropy with parameter a, pA^g) is the re- 
duced density matrix of \G{g))AB with respect to Alice's 
subsystem, an d \\} are the eigenvalues of pA^g) in de- 
creasing order|l6l4l8l|. Thus, if the Renyi entropies of 
two states intercept for some a, they cannot convert to 
each other by LOCC even in the presence of ancillary 
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FIG. 1: Differential local convertibility illustration. Alice and 
Bob control the two bi-partitions of a physical system whose 
ground state depends on a Hamiltonian that can be varied via 
some external parameter g. In one phase (panel a), the con- 
version from one ground state \G{g)) to another \G{g + A)) 
requires a quantum channel between Alice and Bob, i.e, co- 
herent interactions between the two partitions are required. 
We say that this phase has no local convertibility. After phase 
transition (panel b), Alice and Bob are able to convert \G{g)) 
to \G{g-\-A)) via only local operations and classical communi- 
cations, and local convertibility becomes possible. If A ^ 0, 
local convertibility becomes differential local convertibility. 
This implies that it is impossible to completely simulate the 
adiabatic evolution of the ground state with respect to g by 
LOCC in one quantum phase phase and possible in the other 
phase. 



entanglement [19]. In the A ^ 0+ limit, we may instead 
examine the sign of dgSa{g) for all a. U dgSa{g) does 
not change sign, the effect of an infinitesimal increase of 
g results in global shift in Sa{g)^ with no intersection be- 
tween Sa{g + A) and Sa{g)- Otherwise, an intersection 
must exist. 



RESULTS 
Transverse field Ising model. 

Before we consider the general XY model, we highlight 
key ideas on the transverse Ising model, which has the 
Hamiltonian 



N 



(2) 



where cr^, for k = x^y^z are the usual Pauli matrices and 
periodic boundary conditions are assumed. The trans- 
verse Ising model is one of the simplest models that has a 
phase transition, therefore it often serves as a test bed for 
applying new ideas and methods to quantum phase tran- 



sitions. Osterloh et al. have previously shown that the 
derivative of the concurrence is a indicator of the phase 
transition^. Nielsen et al. have also studied concurrence 
between two spins at zero or finite temper at ure[5|. Re- 
cently, the Isin^chain with frustration has been realized 
in experiment j20|. 

The transverse Ising model features two different quan- 
tum phases, separated by a critical point at = 1. When 
^ < 1, the system resides in the ferromagnetic (symmet- 
ric) phase. It is ordered, with nonzero order parameter 
(cr^), that breaks the phase flip symmetry H^crf . When 
^ > 1 , the system resides in the symmetric paramagnetic 
(symmetry broken) phase, such that (a^) = 0. 

There is systematic qualitative difference in the com- 
putational power afforded by perturbation of g within 
these two differing phases. In the paramagnetic phase, 
dgSa{g) is negative for all a, hence increasing the ex- 
ternal magnetic field can be simulated by LOCC. In the 
ferromagnetic dgSa{g) changes signs for certain a. Thus 
the ground states are not locally convertible, and per- 
turbing the magnetic field in either direction results in 
fundamentally non-local quantum effects. 

This result is afforded by the study of how Renyi en- 
tropy within the system behaves. From Eq.(l), we see 
that Renyi entropy contains all knowledge of {A,}(2l|. 
For large a, Xf vanishes when is small, and larger 
eigenvalues dominate. In the limit where a oo, all 
but the largest eigenvalue Ai may be neglected such that 
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log2Ai. In contrast, for small values of a. 



smaller eigenvalues become as important as their larger 
counterparts. In the a ^ 0+ limit, the Renyi entropy 
approaches to the logarithm of the rank for the reduced 
density matrix, i.e, the number of non-zero eigenvalues. 

This observation motivates study of the eigenvalue 
spectrum. In systems of finite size (Fig. 2 (a)), the 
largest eigenvalue monotonically increases while the sec- 
ond monotonically decreases for all g. All the other eigen- 
values Xk exhibit a maximum at some point gk- From the 
scaling analysis (Fig. 2 (b) and (c)), we see that as we 
increase the size of the system, gk ^ 1 for all k. Thus, 
in the thermodynamic limit, A^ exhibits a maximum at 
the critical point of ^ = 1 for all /c > 3. Knowledge of 
this behavior gives intuition to our claim. 

In the ferromagnetic phase, i.e., < g < 1, dgSa{g) 
takes on different signs for different a. When a ^ 0+, 
tends to the logarithm of the effective rank. From Fig. 
2 we see that all but the two largest eigenvalues increase 
with ^, resulting in an increase of effective rank. Thus 
dgSa{g) > for small a. In contrast, when a ^ oo, 
Sa — log2Ai. Since Ai increases with g (Fig. 2), 
dgSa{g) < for large a. Therefore, there is no differen- 
tial local convertibility in the ordered phase. 

In the paramagnetic phase, i.e., ^ > 1, calculation 
yields that dgSa{g) is negative for both limiting cases 
considered above by similar reasoning. However, the in- 
termediate a between these two limits cannot be ana- 
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FIG. 2: The largest 4 eigenvalues and the scaUng analysis, (a) 
The four largest eigenvalues of transverse Ising ground state 
for N=10 case. The red line represents the largest eigen- 
value Ai, and the black line is the second one A2. Note that 
we have artificially magnified A3 (green) and A4 (blue) by 
10 times for the sake of clarity since each subsequent eigen- 
value is approximately one order smaller than its predecessor, 
(b) The scaling behaviors of the maximal points of the third 
eigenvalue. When N ^ 00, the maximum points approach 
to the critical point with certain acceptable error. The black 
dots are the data and the red curves are the exponentially fit 
of these data. The maximum point for the third eigenvalue 
gs = 0.495 X + 1.09177. (c) The maximum point 

for the fourth eigenvalue g4 = —0.082 x exp(— -g-^) + 0.9996. 
The maximums points of smaller eigenvalues have similar be- 
havior. 

lyzed in a simple way. The detail and formal proof of the 
result can be found in the 'Method' section, where it is 
shown that dgSa{g) still remains negative for all a > 0. 
Thus, differential local convertibility exists in this phase. 

These results indicate that at the critical point, there is 
a distinct change in the nature of the ground state. Prior 
to the critical point, a small perturbation of the exter- 
nal magnetic field results in a change of the ground state 
that cannot be implemented without two body quantum 
gates. In contrast, after phase transition, any such per- 
turbation may be simulated completely by LOCC. 

XY model. 

We generalize our analysis to the XY model, with 
Hamiltonian 

H = - ^[^(1 + 7K<Tf+i + 1(1 - 7)<Tf af+i + 5<Tf],(3) 

i 

for different fixed values of 7 > 0. The transverse Ising 
model thus corresponds to the the special case for this 
general class of models, in which 7 = 1. For 7^1, 
there exists additional structure of interest in phase space 
beyond the breaking of phase flip symmetry at ^ = 1. 
In particular, there exists a circle, + 7^ = 1, on 
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FIG. 3: XY model local convertibility phase diagram. Con- 
sideration of differential local convertibility separates the XY 
model in three phases, which we label phase lA, phase IB 
and phase 2. We consider differential local convertibility for 
fixed values of 7 while g is perturbed. Differential local con- 
vertibility is featured within both phase IB and phase 2, but 
not phase lA. 

which the ground state is fully separable. The functional 
form of ground state correlations and entanglement are 
known differ substantially on either side of the circle0- 
01, which motivates the perspective that the circle is 
a boundary between two differing phases. Indeed, such 
a division already exists from the perspective of an en- 
tanglement phase diagram, where different 'phases' are 
characterized by the presence and absence of parallel en- 
tanglement [25]. 

Analysis of local convertibility reveals that from the 
perspective of computational power under adiabatic evo- 
lution, we may indeed divide the system into three sep- 
arate phases (Fig. 3). While the disordered paramag- 
netic phase remains locally convertible, the local con- 
vertibility of the ferromagnetic phase now depends on 
whether + 7^ > 1. In particular, for each fixed 
7 > 0. The system is only locally non-convertible when 
g > Y^l — 7^. We summarize these results in a 'local- 
convertibility phase-diagram', where the ferromagnetic 
region is now divided into components defined by their 
differential local convertibility. 



Different partitions. 

Numerical evidence strongly suggests that our results 
are not limited to a particular choice of bipartition. We 
examine the differential local convertibility when both 
systems of interest is partitioned in numerous other ways, 
where the two parties may share an unequal distribution 
of spins (Fig. 4). The qualitative properties of dgSa re- 
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FIG. 4: The sign distribution of dgS^ in Ising model and 
XY model for different bi-partitions. The sign distribution 
of dgSa on the a — g plane for different bi-partitions on the 
systems. Panel (a), (5), (c) correspond to Ising model of size 
A/" = 12, where Alice possesses 6, 5, 4 of the spins respectively. 
dgSa is negative in lighter regions and positive in the red 
regions. Clearly, regardless of choice of bipartition, dgSa is 
always negative for ^ > 1 and takes on both negative and 
positive values otherwise. Note that for very small dgSa 
only becomes negative for very large a and thus appears com- 
pletely positive in the graph above. The existence of nega- 
tive dgSa can be verified by analysis of dgSa in the a ^ oo 
limit. The choice of bipartition affects only the shape of the 
dgSa = boundary, which is physically unimportant. Panel 
(d), (e), (/) correspond to XY model with fixed 7 = 

= 14 and L = 7, 6, 5 respectively. Panel (/i), (i), (j) cor- 
respond to XY model with fixed 7 = a/7/4, = 18 and 
L = 9, 8, 7, respectively. Here the value of L represents a bi- 
partition in which L qubits are placed in one bipartition and 
N — L qubits in the other. The only region in which dS^/dg 
takes on both negative and positive values is in phase lA of 
Figure 4. Note that the transition between phase lA and 
IB occurs at p = 0.5 for 7 = V3/2 (Point E, in Fig 3) and 
g = 0.75 for 7 = V7/4 (Point D, in Fig 3.) 



main unchanged. While it is impractical to analyze all 2 
possible choices of bipartition, these results motivate the 
conjecture that differential local convertibility is indepen- 
dent of our choice of partitions. Should this be true, it has 
strong implications: The computation power of adiabatic 
evolution in different phases are drastic. In one phase, 
perturbation of the external field can be completely sim- 
ulated by LOCC operations on individual spins, with no 
coherent two-spin interactions. While, in other phases, 
any perturbation in the external field creates coherent in- 
teractions between any chosen bipartition of the system. 



DISCUSSION 

The study of differential local convertibility of the 
ground state gives direct operational significance to phase 
transitions in the context of quantum information pro- 
cessing. For example, adiabatic quantum computation 
(AQC) involves the adiabatic evolution of the ground 
state of some Hamiltonian which features a parameter 
that varies with time [2^, H^l- This is instantly remi- 
niscent of our study, which observes what computational 
processes are required to simulate the adiabatic evolu- 
tion of the ground state under variance of an external 
parameter in different quantum phases. 

Specifically, AQC involves a system with Hamiltonian 
(1 — s)Ho + sHp^ where the ground state of Hq is simple 
to prepare, and the ground state of Hp solves a desired 
computational problem. Computing the solution then in- 
volves a gradual increment of the parameter s. By the 
adiabatic theorem, we arrive at our desired solution pro- 
vided the s is varied slowly enough such that the sys- 
tem remains in its ground state [26, 27]. We can regard 
this process of computation from the perspective of local 
convertibility and phase transitions. Should the system 
lie in a phase where local convertibility exists, the incre- 
ment of s may be simulated by LOCC. Thus AQC cannot 
have any computational advantages over classical com- 
putation. Only in phases where no local convertibility 
exists, can AQC have the potential to surpass classical 
computation. Thus, a quantum phase transition could 
be regarded as an indicator from which AQC becomes 
useful. 

In fact, the spin system studied in this paper is directly 
relevant to a specific AQC algorithm. The problem of 
"2 — SAT on a Ring: Agree and Disagree" features an 
adiabatic evolution involving the Hamiltonian 

_ AT N 

H{s) = {l-s) ^(1 - aj) + . ^ -(1 - (4) 

where s is slowly varied from to 1 This is 

merely a rescaled version of the Ising chain studied here, 
where the phase transition occurs at 5 = |. According 
to the analysis above, the AQC during the paramagnetic 
phase can be simulated by local manipulations or classi- 
cal computations. For the period of ferromagnetic phase, 
we can do nothing to reduce the adiabatic procedure. 

In this paper, we have demonstrated that the compu- 
tational power of adiabatic evolution in the XY model 
is dependent on which quantum phase it resides in. This 
surprising relation suggests different quantum phases 
may not only have different physical properties, but may 
also display different computational properties. This 
hints that not only are the tools of quantum information 
useful as alternative signatures of quantum phase tran- 
sitions, but that the study of quantum phase transitions 
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may also offer additional insight into quantum informa- 
tion processing. This motivates the study of the quantum 
phases within artificial systems that correspond directly 
to well known adiabatic quantum algorithms, which may 
grant additional insight on how adiabatic computation 
relates to the physical properties of system that imple- 
ments the said computation. There is much potential 
insight to be gained in applying the methods of analysis 
presented here to more complex physical systems that 
featuring more complex quantum transitions. 

In addition, differential local convertibility also may 
possess significance beyond information processing. One 
of the proposed indicators of a topological order involves 
coherent interaction between subsystems that scale with 
the size of the system(2i,[2i|. In our picture, such a indi- 
cator could translate to the requirement for non-LOCC 
operations within appropriate chosen bipartite systems. 
Thus, differential local convertibility may serve as an ad- 
ditional tool for the analysis of such order. 



METHODS 
Eigenvalue properties. 

For the transverse field Ising model, the largest eigen- 
value Ai monotonically increases while the second A2 
monotonically decreases for all g (Fig. 2). In the ther- 
modynamic limit all the other eigenvalues increase in 
the g < I region and decrease in the g > I region. 
Moreover, the eigenvalues other than the largest two are 
much smaller than Ai and A2. Therefore we can av- 
erage them when considering their contribution to the 
Renyi entropy. Thus the eigenvalues are assumed to be 
0.5 + 5, 0.5 - e, (e - (5)/(2^ - 2), (e - (5)/(2^ - 2), . . . when 

< 1, and 1 - (5' - e^ e^ (57(2^ - 2), ... , when g > 1, 
where n is the particle number belonging to Alice and 
certainly Bob has the other N — n particles. Because 
of some obvious reasons such as Ai > A2 > A3 • • • and 
all these eigenvalues are positive, and so on, we can de- 
rive the following relations easily: < 6 < e < 0.5,0 < 
05 /dg < de/dg,0 < 5' < e' < 0.5,05' /dg < 0, and 
de' /dg < 0. Then we can prove the main result for each 
phase region. Namely, in the g < 1 phase, dgS^ is pos- 
itive for small a but negative for large a; and for the 
g > 1 phase, dgS^ < for all a. 



and 

dSc. 1 



dg I- a (0.5 + 5)^ + (0.5 - e)^ + (2^ - 2) 



^(0.5 + 6) - - (0.5 - e) - + (^) (- - -)J 
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1 - a (0.5 + + (0.5 - e)" + (2^ - 2) 



x|^[(0.5 + ^)-^ 



\dg 
de 
dg 



K0-5-^r-^-(^r ^li- 



ce) 



In the thermodynamic limit A/" 00, (e— ^)/(2^ — 2) 



When < a < 1, {^^^ ^ ^ 00, 
1 a ^de d5^ 



dSc 
dg 



1 — ae — 5^dg dg 



> 0. 



When a > 1, 



-6 



0, 



~di ^ 1-a^d^^ 0.5 -e^ 



de. 



(7) 



(8) 



As (1 + ^ > 1, and < 05 /dg < de/dg, we 

can see that the solution of dS^/dg = (labeled as a^) 
always exists in the region ^ < ip|a > 1, and dS^/dg 
will be negative as long as a > a^. Moreover we can 
also see that the smaller g is, the smaller 5 and e are, 
and the smaller (l + (e + 5)/ (0.5 — e)) is, and therefore 
the larger ao should be. This explains why we need to 
examine larger value of a to find the crossing when g is 
very small. 

Notice that in the above analysis we used the N ^ 00 
condition in the g < 1 region. We can also find that 
in Fig. 4 for finite A" = 12 there is some small green 
area in the region a < 1 < 1. However, in the above 
analysis of infinite A", this area should be totally red. 
This difference is due to the finite size effect. 



Paramagnetic phase. 



Ferromagnetic phase. 

In the ferromagnetic phase ^ < 1, the eigenvalues are 
0.5 + 5, 0.5 - e, (e - (5)/(2^ - 2), (e - (5)/(2^ - 2), ... . So 
the Renyi entropy 

(5) 



In the paramagnetic phase, g > 1. The eigenvalues are 
1-5' - e^ e^ ^7(2^ - 2), ... . The Renyi entropy 

= log[(l -5'- e'r + {e'r + {r - 2)(^)-].(9) 

And 

dSc, _ _l 1 

dg ~ l-a{l-3' -e')^ + {e')^ + (2^-2){^Y 
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dq\ 



dg ^2^-2^ dg. 
1-6' -eT~^ 



l-k— 



2 1 



-k- 



where d5' / dg.de' / dg < 0; and e7(l - - ^0.^7(1 - 
e' - 6'){2'' - 2) G (0, 1), since Ai > A2 > A3. 

So when a > 1, [(. . . Y'^ - 1] < 0, {...}> 0, a/ (I - 
a) < 0. We have dSa/dg < 0; and when < a < 1, 
[(. . . )"-^ - 1] > 0, {...}< 0, a/(l - a) > 0. We also 
have dSa/dg < 0. Hence, we can obtain that when g > 
1, dSa/dg < for all a > 0. If we consider the full 
condition for the local convertibility which includes the 
generalization of a to negative value, it also can be proved 
easily that the local convertible condition dgSa/oc > for 
all a is satisfied in the g > 1 phase. As for the g < 1 
phase, the sign changing in the positive a already violates 
the local convertible condition, so that we do not need to 
consider the negative a part. In fact, generally speaking, 
for the study of differential local convertibility^ we can 
only focus on the positive a part, because the derivative 
of Renyi entropy over the phase transition parameter will 
necessarily generate a common factor a, which will cancel 
the same a in the denominator. 

To conclude. In the g > 1 phase dSa/dg is negative for 
all a, and in the ^ < 1 p| a < 1 region it is positive, while 
in the ^ < 1 P| a > 1 region it can be either negative or 
positive with the boundary depending on the solution of 
Eq.(4). 
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